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ABSTRACT 



Context. Young populations at Z<Z© are being examined to understand the role of metallicity in the first phases of 
stellar evolution. For the analysis it is necessary to assign mass and age to Pre-Main Sequence (PMS) stars. While 
it is well known that the mass and age determination of PMS stars is strongly affected by the convection treatment, 
extending any calibration to metallicities different from solar one is very artificial, in the absence of any calibrators 
for the convective parameters. For solar abundance, Mixing Lenght Theory models have been calibrated by using the 
results of 2D radiative- hydrodynamical models (MLT-a 2D ), that result to be very similar to those computed with 
non-grey ATLAS9 atmosphere boundary condition and full spectrum of turbolence (FST) convection model both in 
the atmosphere and in the interior (NEMO-FST models). 

Aims. While MLT-a 2D models are not available for lower metallicities, we extend to lower Z the NEMO-FST models, 
in the educated guess that in such a way we are simulating also at smaller Z the results of MLT-a 2D models. 
Methods. We use standard stellar computation techniques, in which the atmospheric boundary conditions are derived 
making use of model atmosphere grids. This allows to take into account the non greyness of the atmosphere, but adds a 
new parameter to the stellar structure uncertainty, namely the efficiency of convection in the atmospheric structure, if 
convection is computed in the atmospheric grid by a model different from the model adopted for the interior integration. 
Results. We present PMS models for low mass stars from 0.1 to 1.5 Mq for metallicities [Fe/H]= -0.5, -1.0 and -2.0. 
The calculations include the most recent interior physics and the latest generation of non-grey atmosphere models. At 
fixed luminosity more metal poor isochrones are hotter than solar ones by A log T e a/A logZ ~ 0.03-0.05 in the range 
in Z from 0.02 to 0.0002 and for ages from 10 5 to 10 7 yr. 
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1. Introduction 

The study of Pre-Main Sequence (PMS) stars is impor- 
. £h ' tant to trace the modalities of star formation in space and 
, time, to date young stellar systems by means of age trac- 
■ ers which do not suffer the uncertainty in the physics of 
d 1 the upper main sequence stars, to derive the initial mass 
function of very low mass stars and brown dwarfs, and to 
understand the modalities of the stellar rotational evolu- 
tion and depletion of light elements. Theoretical tracks and 
isochrones provide an essential tool to understand and in- 
terpret the experimental data currently available for young 
objects. Since most of these sources are located in nearby 
Galactic star forming regions, this research has been gen- 
erally limited so far to solar chemistry or close to it (e.g 
ISiess et ail , [200 0). Now t hat young populat i ons a t Z<Z 
are being examined (e.g. iRomaniello et all |2006| ), it be- 
comes essential to understand the role of metallicity in the 
first phases of stellar evolution, thus demanding further in- 
vestigations at lower metallicities. 

On the other hand, the description of the PMS evolution of 
stars is one of the most delicate tasks in the context of stel- 



lar astrophysics: the location of the theoretical tracks on 
the Hertzprung-Russell (HR) plane depends critically on 
the ingredients used to calc ulate the models, the equation 
of state and opacities (e.g. iD'Antonal . I1993T ) . the bound- 
ary conditions used to match the integration of the in terior 
with the atmospheric structure (Baraff e et al. . 2002t). an d 
the treatment of convectio n dD'Antona fc Mazzitellil . Fl994f ). 



* These evolutionary tracks and 
available in electronic form at 
http : //www.mporzio . astro . it/%7Etsa/ 



isochrones 
a WEB 



arc 
site 



Detailed investigations by iMontalban etldl ( 2004 ) showed 
that, at T c ffS where atmospheric convection is present, its 
modelling plays an important role in determining the radius 
of these structures, and thus the exact excursion of the the- 
oretical track on the HR diagram; this is the most critical 
parameter in the computation down to T c ff~4000K. Below 
this temperature, both molecular opacities and convection 
treatment are the dominant uncertainties. Unfortunately, 
convection is a rather complex phenomenon, that is scarcely 
known from first principles, so that it is commonly treated 
by means of purely local approache s, the most popular of 
which is the Mixing Length Theory (|Bohm-Vitense I . [l958. 
MLT), where all the physical uncertainties are hidden be- 
low the unique parameter a — l/H p , being / the mix- 
ing scale and H p the pressure scale height; among other 
attempts to model convection locally, the Full Spectrum 
of Turbul ence (FST) model bvlCanuto fc Mazzitellil (|l99ll 
CM) and lCanuto. Goldman fc Mazzitellil (|1996l . CGM) has 
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been largely used in recent years. 

A potentially powerful tool to "calibrate" the convective 
model for the description of these evolutionary stages is 
the comparison between the theoretical loci in the HR plane 
and the position of PMS binaries, for which a rather precise 
esti mation of the masses of the two compone nts is possible 
(e.g IStassun et all . 12004 iBoden et all . 120051 ). The results 
obtained so far are rather ambiguous, and evidentiate the 
difficulty to find out a unique de s criptio n of convection, 
holding in all cases. iBaraffe et al.l (2002) found that sev- 
eral masses in PMS binaries demand an MLT parameter 
a = 1.9, but some others, lying in the same region of the 
HR diagram, r equire a much lower efficiency, e. g. a = 1. 
IStassun et al.1 (|2004l ) also found that models with ineffi- 
cient convection in the interior should be preferred. In a 
case, for the binary components of V1174 0ri, they found a 
very g ood agreement with the models by iMontalban et al.l 
(2004), employing the FST convection bot h in the in- 
terior and in the atmospheres computed by iHeiter et al.1 
(2002j). 2D and 3D radiative hydrodyn amical simulations 
should provide mor e reali stic results (jTrampedach et al.l . 
1999; Lu dwig et all I2002D . Based on their 2D computa- 
tions, iLudwig et al.l ([1999) provided a calibration of the 



MLT parameter a for various effective temperatures and 
gravities: the meaning of this calibration is that the param- 
eter describes the average efficiency of convection within the 
whole superadiabatic zone, in the region of T c g and gravity 
explored by the 2D models. This efficiency is not constant, 
but varies with the position of t he star in the HR diagram. 
IMontalban fc D'Antona I (|2006l ) adopted such a calibration 
and computed the corresponding tracks (MLT-a 2D tracks) 
for solar metallicity. Exam of binary location with respect 
to the MLT-a 2ZJ tracks shows the same difficulties of previ- 
ous track sets, although they result in excell ent agreement 
with the components of the binary V773 Tau (jBoden et all 
120071 ). The MLT-a 2D tracks are not in agreement with the 
lithium depletion patterns in young clusters with the so- 
lar system abundance, but this additional feature can be 
attributed to the high metal abundance generally ad opted 
for the solar model (|Montalban fe D'Antona! I2006D . The 
same work shows that the FST non-grey tracks of solar 
composition, with boundary conditions based on model at- 
mospheres in which convection adopts the same FST model 
(jHeiter et al.l . 12002]) , present a striking similarity with the 
results of MLT-a 2 ^ models. 

The scope of this paper i s to ex tend the computations by 
IMontalban fc D'Antona I (|2006| ) to lower metallicities, by 
naively assuming that FST models can approximate the 
results of MLT-a 215 simulations also at lower metallicity. 
As we wish to provide a more extended set of results, while 
the FST model atmospheres are available only at T c g > 
4000K, we adopt a way to extend the computation to lower 
T e g. We then use ML T non grey models based on th e atmo- 
spheric structures by lAllard fc Hauschildt 1 (|1997l . AH97) 
and choose the combination of atmospheric and interior 
convection efficiency that provides results similar to the 
FST tracks at T off > 4000K. 

2. The ATON code 

The evolutionary sequences presented in the following 
sections were calculated by means of the ATON code for 
stellar evolution; a detailed descri ption of the numeric al 
structure of ATON can be found in lVentura et all (|2007ft . 



The micro-physics adopted was recently updated, for what 
concerns the radiative and conductive opacities, and the 
equation of state (EOS). 

At low temperatu res, we use th e lates t release of the 
opacity tables by iFerguson et al.1 (|2005l ). completed for 
T > 10000K by the OPAL opacitie s, in the version 
documented by llglesias fc Rogers! (fl99H . 
The EOS adopted in most of the T — P plane is the latest 
release of the OPAL EOS3, overwritten in the pressure 
ionization regime by the EOS from iSaumon et aLT ([1995), 
and extended to the high-density, hig h temperature do- 
main according to the treatment by IStolzmann fc Blockerl 

(pool . 

The formal borders of the convective regions were found 
by means of the classic Schwartzschild criterium. The 
FST scheme (CM or CGM) was used to determine the 
temperature gradient in zones unstable to convection. 
Nuclear burning within convective regions was treated 
according to the instantaneous mixing approximation. 
We computed models for three metallicities [Fe/H] = 
-0.5,-1.0 and 2.0 with an adopted helium mass fraction 
Y=0.25. The solar metallicity is assumed to be Z=0.02 
and the solar mixture f or op acities and EOS is taken 
from iGrevesse fc Sauva l (1999). In this work we adopted 
solar-scaled mixtures in all cases, though a— enhanced 
mixtures are also available. 



2.1. Atmospheric structure and boundary conditions 

In the atmosphere, the opacities have a strong dependence 
on the frequency, so that an atmospheric integration based 
on a Rosseland-type average or on an approximate T(t) 
relation, like in a grey atmosphere, is not adequate for the 
description of the structure. Often this problem is attacked 
by adopting the stratification temperature vs. pressure by 
integrating appropriate model atm ospheres , as fo r exam- 
ple by means of the Kurucz code dKuruczl . Il998f) . Below 
~4000K, the role of triatomic molecules becomes impor- 
tant, and the most adequate model a tmosp heres are so far 
the models by lAllard fc Hauschildtl (|19971 ). Generally, the 
boundary conditions for the interior structure computation 
are then the physical quantities deriving from such integra- 
tion down to a fixed value of the optical depth r, namely 
T p h- Tables of boundary conditions at the chosen r p h are 
used to derive the stellar gravity and T e ff by interpolation. 
This procedure hides a problem: the stellar convection in 
the atmosphere is computed by assuming an efficiency of 
convection. In the MLT model, for example, the grids of 
model atmospheres are computed by fixing the a parame- 
ter to a value a atm . For the interior computation, it may be 
necessary to adopt a different a. For example, in order to 
fit the solar model (that is, to obtain t he solar radiu s at th e 
solar age), a value a; n t = 1.9 is used by IBaraffe et al.l (|2002l ). 
while the AH97 grid adopted for the atmospheric integra- 
tion has a a tm=l- As the MLT must be generally understood 
as a way of obtaining an "average" efficiency of convection, 
more than a model that allows to derive the correct temper- 
ature atmospheric structure, this problem should not worry 
too much, but we should remember that in this way we are 
introducing a dependence of the structure on another pa- 
rameter, namely r p h. This problem is fully discussed by 
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Fig. 1. Evolutionary tracks for the labelled metallicities and for masses 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 
1.1, 1.2, 1.3, 1.4, 1.5 Mq. The solid ones are calculated using NEMO grids of atmosphere models and FST treatment 
of convection in the interior, while dashed lines are MLT-tracks aj, lt =2.0, r p h=3) obtained using AH97 models of atmo- 
spheres. To facilitate the reading of this figure in each panel the evolutionary track corresponding to 0.6 M Q is labelled; 
for this value of mass both the MLT-AH97 and FST-NEMO tracks are reported (see also Fig. H}. 



Montalban et alj (l2004f) . 



Heiter et al.l (|2002f ) made available grids computed by 



means of an improved ver sion of Kurucz code (NEMO; 
iKuruczl 119981: ICastellil . I1997D . They considered both MLT 
models with a atm =0.5, FST models following CM and FST 
models according to CGM. If FST convection is adopted 
both in the atmospheric grid and in the interior, the model 



computa tion does not show temp erature gradients discon- 
tinuities (|Montalban et al l . [2004] ) . 

Different sets of boundary conditions were considered in 
this work: FST-NEMO grids (with CGM convection) and 
the AH97 grids. The boundary conditions (BCs) at fixed 
values of photospheric optical depth r p h contain, for each 
T e ff, gravity and metallicity, the temperature, pressure and 
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(120041), which were pe rformed using the NEMO-FST grids 
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Fig. 2. Non-grey evolutionary tracks for 0.6 M Q and 
Z=0.002 (upper panel) and Z=0.0002 (bottom panel) . The 
solid line is FST tracks with NEMO-FST atmospheres with 
r p h=10. The others are MLT tracks obtained using the 
atmosphere models by AH97 with different values of a m t 
in the computation of the sub-atmospheric convection and 
with the match point at r p h=3. 



geometrical depth at the chosen r p h- The latter quantity is 
the difference between the radius where convection begins 
and V ra d becomes larger than V <2 and the radius at T p h, a 
quantity necessary for the computation of the FST fluxes. 
As suggested by Heiter et al. (2002), we choose r p h=10, to 
avoid discrepancies in the physical quantities due to the tur- 
bulent pressure (not included in the atmosphere modeling 
but present in the integration of the interior) and to reduce 
the differences due to different opacity tables used at both 
sides of r p h. The boundary condition for the internal struc- 
ture are determined by spline interpolation of these tables. 
From initial T e ff and log L/Lq we determine P and T at 
the last point of the internal structure (T p h). An iterative 
procedure is performed until the P(r p h) and T(r p h) values 
derived from the interior and from the atmosphere models 
converge. The T g range of these grids is between 4000- 
10000 K, at lower T e fi we have used AH97 models, which 
include the contribution of many more molecular lines dom- 
inating the opacities with the respect to ATLAS9. In the 
AH97 models, convection is treated with MLT and a a t m =l. 
For metallicities lower than sola,r the available models have 
3000<T cff < 10000 K and surface gravity from log#=3.5 to 
6.0. 



3. The evolutionary tracks 

In Fig. [1] we present evolutionary tracks computed with non 
grey atmospheres for masses in the range 0.1 < M/Mq < 
1.5, and for metallicities [Fe/H]= 0.0, -0.5, -1.0, -2.0. 
In the top-left panel the solar tracks at M > 0.6 Mq are di- 
rectly taken by previous computation by iMontalban et aD 



bv lHeiter et~aT1 (|2002l ) available only for T cff >4000 K. The 
comparison of these solar metallicities tracks with other 
models present in the literature, in particular with those of 
Siess et al. 2000, was made by Montalba n et al. 2004 . 
However IMontalban &: D'Antonal {2006) have shown that 
these models are very similar to the MLT-a 2D models; here 
we extend the computation to lower metallicities. 
Since measurements suggest a primordial abundance of 
deuterium ~ 3.0-5.0 ■ 10" 5 (e.g. lTosil . [l996h . we fix the ini- 
tial deuterium abundance at Xd=4 ■ 10~ 5 by mass, fairly 
representative of the D-abundance of the Population II 
stars. A value ■ 10~ 5 is of t en ad opted for models 

with solar abundances (|Siess et al.l . |2000[) . 
NEMO models are available only for T eff >4000K. Within 
the context of non grey modelling, the only atmospheres 
availa ble at lower T e g are those by lAllard k, Hauschildt I 
(|1997[ ). that use the MLT treatment of convection. We 
therefore decided to extend our models to lower masses 
by using these atmosphere's models. The same convec- 
tion model is adopted in the interior computation, and the 
match with atmospheres is done at r p h=10 in the first case 
and T p h=3 in the second one. The free parameter a- ln t was 
set in order to provide a reasonable continuity between the 
FST tracks and the MLT ones at temperatures just ex- 
ceeding 4000K. Fig [2] shows that ai„t = 2 is a reasonable 
choice especially at lower metallicity (see also the evolution- 
ary tracks for M=0.6 M in FigQJ. For these MLT mod- 
els we choose r p h=3 which is a good co mpromise since we 
avoid having a large influence from a a t m ([Montalban et all , 
2004). The evolutionary tracks computed with the MLT 
treatment of convection are indicated with dashed lines in 
FigUl since the AH97 atmospheres are available only in the 
range of gravities log g > 3.5, we had to skip the deuterium 
burning phase from our computations. 
As outlined in the introduction, PMS stars can be very 
useful age tracers. Comparing the theoretical results with 
the observed loci of stellar sources in associations requires 
the computation of the isochrones. In Fig|3]we report three 
different groups of isochrones (corresponding to 10 5 , 10 6 , 
10 7 yr) for each of the metallicities investigated. 
The youngest isochrones do not include the low mass stars, 
because, as already stressed, we started the evolution after 
the deuterium burning phase. As for very low mass mod- 
els (those with AH97 atmospheres) we skip the D burning 
phase, the ages of these PMS stars must be taken with cau- 
tion. At 10 7 yr, for example, including D burning would have 
increased the ag e by about only 3% at 0.5 M Pl} but by about 
30% at 0.1M© (jD'Antona fc Mazzitellil . 119971 ). Young ages 
as low as 10 5 yr are also uncertain, because of the unsure role 
played by the protostellar accretion phase (jPalla I . I2001I ). It 
is evident from the figure that at fixed age and luminosity 
more metal poor isochrones are hotter by an amount which 
is almost independent of luminosity and that these amounts 
became lower at lower metallicity. The effect is mostly due 
to the opacity reduction, as it is well known in the study of 
main sequence models. 

4. Conclusions 

A grid of stellar evolutionary tracks for low metallicity 
Pre-Main Sequence stars with masses between 0.1 and 1.5 
M Q was presented. These models are based on up-to-date 
physics and updated non grey atmosphere models were 



Di Criscienzo et al.: Updated PMS tracks at low Z for 0.1 < M/M Q < 1.5 



5 




4 3.9 3.8 3.7 3.6 3.5 



Log(T e „) 

Fig. 3. Comparison between the isochrones calculated for 
10 5 , 10 6 , 10 7 yr for different metallicities (from left to right: 
[Fe/H]=-2.0, -1.0, -0.5, 0) 
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used. A coherent treatment of convection in the interior and 
exterior region of the star was employed at T ff >4000K. We 
extended our computations to models of smaller masses by 
using the AH97 grid of model atmospheres. The parame- 
ters T p h and ai n t were chosen to provide a smooth transition 
between the two model sets. This, of course, is only an ed- 
ucated guess to the problem of Pre-Main Sequence models 
at low metallicity. 

The models (available in electronic form at the WEB loca- 
tion http://www.mporzio.astro.it/%7Etsa/) can now 
be confronted to the complex realm of very young objects, 
providing important information on ages and star forma- 
tion processes, and, on the other hand, providing some new 
constraints for PMS models. 
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